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Abstract 

The gas of interacting excitons in quantum wells is studied. We obtain 
the Hamiltonian of this gas by the projection of the electron-hole plasma 
Hamiltonian to exciton states and an expansion in a small density. Matrix 
elements of the exciton Hamiltonian are rather sensitive to the geometry of 
the heterostructure. The mean field approximation of the exciton Hamiltonian 
gives the blue shift and spin splitting of the exciton luminescence lines. We 
also write down the Boltzmann equation for excitons and estimate the energy 
and spin relaxation time resulting from the exciton-exciton scattering. Making 
use of these calculations we succeeded to explain some recent experimental 
results which have not been explained so far. 
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I. INTRODUCTION 



Optical properties of dense exciton gas in quantum wells had been studied intensively 
both experimentally and theoretically in the past decade. Originally, attempts to obtain 
dense exciton gas were motivated by a desire to reach Bose - Einstein condensation of 
excitons. However, the fascinating observed physical phenomena and the promising potential 
for device application had made the exciton gas a very interesting system by itself. 

Recent experiments have stressed the importance of exciton-exciton interaction of two- 
dimensional (2D) excitons in quantum wells. Exciton-exciton interaction affects exciton 
photoluminescence, breaks the symmetry between excitons with different spin components, 
contributes to depolarization of the exciton gas, and controls spin and momentum relaxation 
in the exciton system. A thorough investigation of exciton-exciton interaction is crucial for 
understanding of a number of experimental results in exciton gas. 

Blue shift of the exciton luminescence line with growing exciton density was reported in 
single quantum wells , multiple quantum wells 0,|4| and coupled quantum wells where 
electrons and holes are spatially separated The blue shift is attributed to the net 

exciton-exciton repulsion interaction in quantum wells J?| [|. 

Other interesting interaction induced phenomena are related to exciton spin. Here we 
use the term spin for the projection of the total angular momentum of the exciton to the 
direction perpendicular to quantum well plane (the growth direction). In quantum wells 
size quantization leads to an energy separation between heavy and light holes. Therefore 
one can treat the ground state exciton as a bound state of a conduction electron with spin 
±1/2 and a heavy hole with spin ±3/2. So, in the ground state, i.e., s-state exciton, the 
spin can take the values ±1 and ±2. Total momentum conservation allows only states with 
spin projection ±1 to be optically active in single photon experiments. 

Spin ±1 and spin -1 excitons can be created independently by pumping circular polarized 
light. A spontaneous energy splitting between the two different components of spin-polarized 
exciton gas in zero magnetic field was observed in multiple quantum wells [|T0|-|T3"| and in 
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coupled quantum wells |L4j]. The splitting increased with the exciton density and decreased 
with the separation of the electron and hole induced by external electric field in coupled 
quantum wells. The splitting was explained by the exciton-exciton exchange interactions 

At high exciton densities the decay of the exciton luminescence cannot be characterized 
by one decay time. Initially the decay is very fast and then it is followed by a slower 
relaxation. |Jl5Hl7| Similar results were obtained for the decay of luminescence polarization. 
T6| , [T7|| A very careful analysis of the luminescence intensity and polarization by Baylac, 



Amand et al. ]16]|T7| and two-pulse experiments performed by Le Jeune et at |T2| proved 
that the short decay time can be explained only by exciton-exciton scattering. In dense 
exciton gas exciton-exciton scattering can also be the leading mechanism of momentum 



relaxation. 18 



The active investigation of exciton-exciton interaction and many new experimental re- 
sults makes it difficult to overestimate the importance of theoretical description of this 
interaction. Such a description, however, encounters significant difficulties even when the 
exciton density is small and the interaction between excitons could be considered with the 
help of perturbation theory. Indeed, let us consider two electron - hole pairs (el, hi) and 
(e2,h2) bound in excitons. One could expect that the Coulomb interaction within pairs 
(i.e., the interaction between el and hi and the interaction between e2 and h2) has to be 
taken into account exactly because it provides bound states, while the Coulomb interaction 
between particles belonging to different excitons (e.g., the interaction between el and h2 
and the interaction between e2 and hi) can be considered as a perturbation. However, 
due to electron-electron and hole-hole exchange it is impossible to say if bound pairs are 
really (el,hl) and (e2,h2) or (el,h2) and (e2,hl). For this reason it is not clear what part 
of the Coulomb interaction is the interaction between excitons and can be considered as a 
perturbation. 

A few approaches have been developed which formally related the gas of electron-hole 
pairs to a Bose gas and avoided such difficulty. Hanamura [19| made use of the Usui trans- 



formation, |20| that makes a correspondence between the space of fermion pairs and a "hy- 



pothetical" boson space. The transformed Hamiltonian still contained Fermi operators. To 
eliminate them the commutators of this Hamiltonian with Bose operators averaged over the 
ground state of the fermion system were declared to be equal to the commutators of the 
target boson Hamiltonian with Bose operators. The harmonic part of the resulting boson 
Hamiltonian immediately led to the known exciton states. Thus, it was natural to assume 



that the anharmonic part describes exciton-exciton interaction. Haug and Schmitt-Rink |21 
(see also Ref. |§) modified this approach introducing creation and annihilation operators as 
linear combinations of pairs of Fermi operators. Commutation relations of new operators 
are different from those of Bose operators by terms proportional to the density of bosons. 
The coefficients of the linear combination were found from a variation principle and ap- 
peared to be single-exciton wave functions corrected by interactions between bosons. Stolz 
et al. |22| used another approach where they represented wave functions of the electron-hole 



plasma as linear combinations of products of single-exciton wave functions. Then again a 
variational principle was used to find the single-exciton wave function. In all these works 
the exciton-exciton interaction was calculated in the leading order in the exciton density. 

In this paper we develop a theory of exciton-exciton interaction in quantum wells which 
is free of a formal definition of exciton creation and annihilation operators. We do not use 
a variational method, the accuracy of which is difficult to control. The physical basis to our 
approach is that we consider a system of equal number of electrons and holes at low enough 
temperature and low density, when all the particles are bound in excitons. If the density is 
not small then excitons overlap and the Coulomb interaction between particles of different 
excitons is of the same order as the electron-hole interaction within one exciton. In such a 
case the identity of excitons is lost and the electron-hole plasma can hardly be considered as 
gas of excitons. If the other limitation, i.e., small kinetic energy of excitons is not met then 
excitons can be ionized and the system can not be described as the gas of excitons only. 
Experimentally small exciton kinetic energy is achieved by resonant excitation at lattice 
temperature much lower than the binding energy []B|, |T~2| , [T5| p77Tl . 



Technically our approach resembles that of Stolz et al. [p2jl - However, we don't rely on 
a variational principle but use a systematic expansion in small exciton density. Finally, we 
arrive at an expression for the second quantized exciton Hamiltonian with coefficients that 
are expressed in terms of the free single-exciton wave function. This function, and so the 
coefficients, strongly depend on the geometry of the heterostructure. 

We use the exciton Hamiltonian for the calculation of the blue shift of the exciton lumi- 
nescence line and the energy splitting between excitons with different spin projection. These 
phenomena are described by Hamiltonian matrix elements which are diagonal with respect 
to occupation numbers of exciton states. Off-diagonal matrix elements of the Hamiltonian 
describe exciton-exciton scattering. We use them to write down the Bolzmann equation for 
the excitons and estimate the relaxation time in the exciton gas. For numerical calculation 
of the coefficients of the Hamiltonian we use the variational single-exciton wave function 



that we suggested earlier [23 



We are not trying to reach an exact numerical matching of our theory with experiments. 
This would require an accurate calculation of the single-exciton wave function for different 
geometries. Our goal is rather to explain semi-quantitatively (i.e., with an accuracy better 
than the order of magnitude) as many experiments as possible. In other words, we are trying 
to show that our theory is able to describe all so far detected phenomena in the exciton gas 
of a small density which are related to exciton-exciton scattering. 

The structure of the paper reflects the specifics of the problem. We consider electrons 
and holes confined in quantum wells. Their motion in the growth direction is strongly 
quantized (we assume that only the ground electron and hole states are occupied) and 
the problem is essentially two-dimensional. So, in the next section we reduce the three- 
dimensional (3D) Hamiltonian of electrons and holes to a 2D one. In Sec. [TTI| we present 
the main assumptions of our approach and give the derivation of the exciton Hamiltonian. 
Some cumbersome details of the derivation are transferred to appendices. In Sec. [TV] we 
calculate the blue shift of the exciton luminescence line and the energy splitting of excitons 
with opposite spins. In Sec. |V] we estimate the exciton-exciton relaxation time with the 



help of the Bolzmann equation. We discuss our results and compare them to experiments 
in Sec. [VI]. Our conclusions are presented in Sec. |V11| . 



II. THE HAMILTON! AN OF 2D ELECTRON - HOLE PLASMA 

Due to a strong quantization of the electron and hole motion in the growth direction the 
dynamics of electrons and holes is essentially two-dimensional. For this reason excitons in 
quantum wells sometimes are considered as purely 2D [|7|-|J . In this approximation, however, 
it is impossible to describe effects of the geometrical parameters of quantum wells on the 
exciton binding energy and their interaction. So we use a more realistic model which takes 
into account the geometry of the heterostructure. 

In quantum wells, light and heavy holes are split off in energy. We consider only heavy 
holes, assuming that the splitting is much larger than kinetic energies of all involved particles 
and the interaction between them. Therefore, light hole states are not occupied and the 
Hamiltonian of the plasma of N electrons and N heavy holes in quantum wells is 

= E i T ej + T hj + U 0e (z ej ) + U 0h (z hj )] (2.1) 
3 

+ Y U eh (\f ei - f hj \,Z ei , Z h j) + - Y Pee(\ ) + Uhh{\^hi ~ r hj\i z hh z hj)] ■ 

ij ij 

Here r e and z e are electron coordinates in the quantum well plane and in the growth direction 
respectively, r h and Zh is the same for holes, 

h 2 v 2 h 2 d 2 

T e = - — ^r~o , (2.2a) 

2m e 2m e ozj 

h 2 X7 2 f) 2 f) 2 
Th = -*JLh - _!!_^_ , (2.2b) 
2m\\ 2m± oz^ 

are the operators of the electron and hole kinetic energy, m e is the electron mass, m\\ and 
m± are in-plane and perpendicular hole effective masses, V e and are the derivatives with 
respect to in-plane coordinates of electron and hole, U 0e (z e ) and U h(zh) are the heterostruc- 
ture potentials that confine electrons and holes in quantum wells, and U e h, U ee , and Uhh are 
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the electron-hole, electron-electron, and hole-hole Coulomb interaction energies respectively. 
If the difference between dielectric constants in wells and barriers is negligible then 

e 2 

Uehdn - r 2 \, zx, z 2 ) = = = , (2.3a) 

K\J (Zi - Z 2 ) + (ft - f 2 ) 

e 2 

Ueedn -f 2 \,z 1 ,z 2 ) = Uhhdn -f 2 \,z u z 2 ) = — = = , (2.3b) 

KyJ (zi - z 2 ) + (f x - f 2 ) 

where k is the dielectric constant. 

The method that we develop can be applied to a single quantum well or to coupled 
quantum wells with all possible well and barrier widths. In this paper we use it only for 
coupled quantum wells where electrons and holes are spatially separated. We neglect the 
electron-hole exchange p4| , p5| which is very small in general [^,[27]] and is further reduced 



in coupled quantum wells due to the smaller electron-hole wave function overlap. We also 
neglect the reduced symmetry of quantum well - barrier interfaces since this effect is very 
small. |28| As a result, Hamiltonian ( |2.1| ) does not depend on electron and hole spins. 



The Hamiltonian ( |2.1| ) contains an information about particle confinement in quantum 
wells which makes the study of in-plane states of the particles more difficult. Fortunately, in 
the most interesting cases the in-plane motion can be separated and the problem is simplified. 
Typically the interaction energy between the lowest level in a quantum well and the first 
excited one (a few tens of meV or even larger) is much larger than Coulomb interaction 
between particles (a few meV). For this reason at low enough temperature particles are 
confined at the lowest quantization level. A distortion of the wave functions caused by an 
admixture with the next level due to Coulomb interaction between the particles can be 
neglected. In such dependence of electron and hole wave functions is described by 

Ce( z ) and Ch( z ) respectively, which satisfy the equations 

h 2 d 2 C 

"2^T + Uoe{z)Ce = EoeCe ' (2,4a) 

+ U oh {z)( h = E oh ( h . (2.4b) 



h 2 d 2 ( h 



2m_i_ dz 2 

Here E Qe and E h are the electron and hole confinement energies respectively. The wave 



function of the gas is a product of all single particles wave functions ( e and (h arid a many- 
particle wave function describing the in-plane state. The Hamiltonian that controls the 
in-plane dynamics of the gas is obtained from the Hamiltonian, ( |2.1|) , by averaging with the 
functions ( e and (h, and has the form 



n 2 v% n 2 v 2 h ; 



2m P 2mw 



+ ^2 u eh(\r e i - r hj \) + - [u ee (\f ei - f ej \) +u hh (\f hi - f hj \)] , (2.5) 



where 



U; 



v( r ) = / Uij(r, z 1} z 2 )( i (z 1 )( j (z 2 )dz 1 dz2 . (2.6) 



All details of the derivation of 2D Hamiltonian, ( |2.5| ), from 3D one, ( |2.1| ), in the case of two 
particles are given in Ref. [ 23| . The derivation of many particle Hamiltonian can be done in 



the same way. 

Thus to study electron - hole gas confined in quantum wells it is enough to consider only 
a 2D problem described by Hamiltonian, (|2.5|) . 



III. DERIVATION OF THE EXCITON HAMILTONIAN 

In this section we derive the 2D Hamiltonian of exciton gas starting from the 2D Hamil- 
tonian of electrons and holes, Eq. ( |2.5| ). We begin this derivation with a discussion of the 
necessary conditions for considering the electron - hole gas as exciton gas. 

Not in any state the electron - hole gas can be represented as a gas of excitons. To make 
this possible, two conditions are necessary. The first one is a small enough concentration of 
electrons and holes, n. Under this condition all electrons and holes can be bound in excitons 
and the excitons are far from each other. That is 

no 2 < 1 , (3.1) 

where a is the exciton radius. If this condition is not met then excitons overlap and the 
Coulomb interaction between electrons and holes of different excitons becomes of the order 
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of the interaction within one exciton. In such a case it is impossible to identify excitons, 
and the electron - hole plasma hardly can be described as a gas of excitons. 

The other condition is that typical exciton kinetic energy is much smaller than the 
absolute value of the exciton binding energy, e&, 

h 2 K 2 



< e b . (3.2) 



2M 

— * 

Here K and M are the exciton wave vector and the exciton mass respectively. If this 
condition is not met then as a result of collisions between excitons they can be excited from 
the ground state or even be ionized into free electrons and holes. Then it is necessary to 
consider excitons interacting with each other and with free electrons and holes. [2Jj We are 



not going to consider this case. Strictly speaking, the accurate condition for ground state 
excitons contains the energy separation between the ground state and the first excited state 
of exciton in the right hand side of Eq. ( |3.2| ). But this difference is of the order of e& and 
we can use Eq. ( |3.2| ) as it is. 

To derive the exciton Hamiltonian we construct a basis for the space of all electron 
- hole states from products of single exciton wave functions. Under the conditions ( ft.lp 
and ( |3.2|) we expect that only the ground state exciton wave function is important. So of 
all matrix elements of Hamiltonian (|2.5|) only the elements between the states constructed 
of these functions should be kept. Such a program in an accurate form is carried out in 
Sec. |III A| , where we show how to calculate all matrix elements of the exciton Hamiltonian. 



Due to inequality fl3.1|) some of these matrix elements are small. Namely, matrix elements 
describing triple and higher order interaction between excitons are proportional to higher 
power of the exciton concentration than matrix elements describing pair interaction, and 
we neglect them. For the calculation of the pair interaction matrix elements it is enough 
to consider the system of only two excitons. This calculation is performed in Sec. |111 B 



Finally, in Sec. [Ill (J| , we use these matrix elements to write down the Hamiltonian of the 
exciton gas in the second quantized form. 



9 



A. Reduction of the electron - hole Hamiltonian to the exciton Hamiltonian 



1. Exciton basis 

The starting point of our derivation is the Hamiltonian of many electrons and holes, 
( |2.5| ). The general idea of the derivation is that we reduce the space of all electron-hole 
states to the subspace of ground state excitons states only, and project the Hamiltonian to 
this subspace. The first step in this direction is the construction of a basis in which this 
subspace can be separated. 

We construct such a basis of symmetrized products of single exciton wave functions. 
Single exciton wave functions are eigenf unctions of the Hamiltonian 

U leh = -— ^ - — ^ + u eh {\r e - r h \) . (3.3) 
2m P 2mu 



Since the Hamiltonian (|3.3|) is independent on electron and hole spins, a single exciton wave 
function can be written as 

^Kas^e, <T e \ f h , tTfc) = g s (<r)l/jg a (f e \ f h ) ■ (3.4) 



Here, ipgai^e', ^h) is the eigenfunction of Hamiltonian, (|3.3|) , and the spin wave function 



g s (a) = 6 SjCT , (3.5) 

can be represented as the product of the electron and hole spin functions, g s (&) = 
g Se (a e )g Sh ((Th)- The projection of electron and hole spins to z-direction can take the values 
s e = ±1/2 and Sh = ±3/2 respectively. The exciton spin projection to the same direction 
is s = s e ± Sh and it can take the values ±1, ±2. The spin variables, a = a e ± cr^, take the 
same values. There is one to one correspondence between the set of electron and hole spins 
and the exciton spin, i.e., each exciton spin corresponds to a single combination of electron 
and hole spins. 

The in-plane exciton wave function is 

lfo a fo r h ) = -Le^ Q (|f e - 4|) . (3.6) 
10 



Here, R = (m e r e + muff^/M is the exciton center of mass coordinate, M = m e + mu is the 
exciton mass, S is the normalization area, and the function, Q (r) is an eigenfunction of the 
Hamiltonian 

h 2 v 2 

Hi eh = + u eh (r) , (3.7) 

where 

+ (3.8) 
/i m e m\\ 

is the reduced mass. 

The functions ^ g a s (r e ,a e ;rh, Ph) form a complete basis for electron-hole pair states. 
That means that a complete basis for the gas consisting of N electrons and N holes can 
be constructed of these functions. This basis consists of correctly symmetrized products of 



single exciton functions, [22 



1 N 

®{ u }(r e u °"ei; r h i, a hl ; . . . , f eN , a eN \ r hN , a hN ) = — ^(-l) p J| ^ Uj (r ejl , a ejl ; f hj2 , a hj2 ) . 

iV ' 3=1 

(3.9) 

Here v stands for the set of quantum numbers (K, a, s), {u} = u%, u 2 , . . . , vjy, the summation 
is carried out over all transpositions of electrons and holes, j± and and P is the parity 
of a transposition. Basis (|3.9| ) is very convenient for our purpose because it easily allows 
us to separate the subspace of wave functions containing only ground state excitons. This 
subspace contains those &{ u } in which all a correspond to the ground state. We enumerate 
this subspace as 1 and the subspace of functions containing at least one excited exciton as 
2. 

It is necessary to note that the functions &{ u } are not orthogonal, in spite of the orthogo- 



nality of the single exciton functions ty v . |?z] Let us, for instance, consider the integral of the 
product of two functions and Q{ v '} which differ by only one of all the quantum numbers 
v, e.g., V\ in the first function is replaced by u[ ^ v\ in the second one. Comparing these two 



wave functions we see that for each term in the sum (|3.9| ) of one wave function there are terms 
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in the sum of the other that differ by exchange of pairs of electrons and holes in such a way 
that the identity of excitons is not violated [e.g., ^ V1 (r e i, cr e i', rki, c^i)^^^^, <7 e2 ; fh 2 , cr^ 2 ) 
and ^/y/ (f e 2, o e2 \ rt2, cr/i2)^i/ 2 (^*ei, cr e i] ^hi, o'hi)]- The integrals of the products of these terms 
equal zero because of the orthogonality of single exciton wave functions. But there are also 
terms that differ by a transposition of electrons or holes which violates the exciton identity, 
e.g., *^(r e i, a e i] r h i, (T hl )^f U2 (f e2 , a e2 ; f h2 , o h2 ) and (r el , a el ; r h2 , ^ 2 )*„ 2 (r e2 , a e2 ; r hl , a hl ). 
The integral of the product of these terms is nonzero. From the definition of the single 
exciton wave functions, Eq. (|3.6|) , we see that this integral, contains the factor a 2 /S. 

If functions arid differ by more than just one value of v or the transposition 
mixes more than a pair of excitons, the integral of contains the factor of a 2 / S to 

a higher power. Eventually, in the calculation of observable quantities each factor of 1/S is 
accompanied by a sum over occupied states of the system which gives the factor of N. So 
that the nonorthogonality of the basis is characterized by the parameter na 2 . 

The basis (|3.9| ) is not normalized. The normalization integral of any of the basis functions, 
f \&{u}\ 2 T\j d 2r ejd 2r hj, contains integrals of the same type of exchange that leads to the 
nonorthogonality. Thus this integral is different from unity by terms of the order of na 2 . 

2. Wave equation for excitons 

After the characterization of the basis (|3.9|) we begin the derivation of the exciton Hamil- 
tonian. We write down the Schrodinger equation for the gas of N electrons and N holes, 
7m> = Ety, in the matrix form in basis ( |3.9|) . To do this we need to represent the eigenfunc- 
tion \1/ as an expansion in the functions and to obtain equations for the coefficients 

of this expansion. For this purpose we multiply the Schrodinger equation by the functions 
complex conjugated to and integrate over all variables. To write down the result in 
the matrix form we introduce the notation for the column of the expansion coefficients 
of the basis function belonging to subspace 1 (only ground state excitons) and the notation 
\l/ 2 for the column of the expansion coefficients of the basis function belonging to subspace 
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2 (including also excited states). Then the matrix equation takes the form 

{H u - N 11 E)m l + (H 12 - N l2 E)m 2 = , (3.10a) 
(H21 ~ N 2 iE)m l + {H 22 - N 22 E)^ 2 = , (3.10b) 

where Tiij are the matrices with elements (${*,}|7i| < l ) { ! /}) and M%j are matrices with elements 
where and belong to i and j subspaces respectively. We use Eq. 

( p.lObQ in order to express ^/ 2 in substitute the result in Eq. (|3.10a|) and come up with 
the equation 

(Hn - Mi5)fi - H iexex) ^i = , (3.11) 
where the effect of excited states is described by 

n (e Xex ) = (n w -M 12 E)('Hz2-M22E)- 1 (H2 1 -N'2iE) . (3.12) 



Although Eq. ( 3.1 1|) contains only it does not have the form of the Schrodinger equation 



because 7-[( exex ) is a nonlinear function of the energy E. The reason is that Eq. ( 3.11 ) is 



equivalent to Eq. ( |3.10 ) and describes the general situation, where electrons and holes can 



occupy excited states. To describe a system where all electrons and holes are bound in 
ground state excitons we make use of the small parameters ( |3.1|) and ( |3.2|) . 



3. Exciton Hamiltonian 



Here we start from Eqs. Q3.ll ) and ( 3.12 ), neglect small terms and obtain the Hamilto- 



nian and the Schrodinger equation for the exciton gas. 

First we note that Afn and Tin contain terms of different order in a 2 /S. The es- 
timate of both Aii and Tin can be done in the same way. The larger the overlap 
between &{ u } and ${1/}, the larger the matrix element (&{ u } is- As we ex- 

plained above, matrix elements between states that differ by a transposition of two 
electrons or two holes contain the factor a 2 /S and describe two-exciton interaction in 
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Tin- Matrix elements between the states that differ by a transposition of three elec- 
trons or three holes, e.g. V ui {r el , a el ; f hl , cr h i)^ U2 (f e2 , cr e2 ; r h2 , a h2 )^ V3 {r e ^ cr e3 ; r^a^) and 
#^(r e i,cr e i;f ft2 ,<XM)#^(r e 2,^ contain (a 2 /S) 2 . \nHu they 

describe a triple exciton interaction that is not reduced to pair interaction. We neglect such 
and other high order terms both in A/n and Tin- 

Matrix elements of A/12, and A/21 are nonzero only because of the nonorthogonality of 
basis functions (|3.9|) . As we showed above (&{ u }\§{ u ty) ~ a 2 / S when v 7^ v' . The matrix 
elements of 7i\ 2 and 7^21 contain the same parameter. To prove this, let us examine the 
largest matrix elements of 7ii 2 . The maximal overlap between subspace 1 and subspace 2 is 
achieved when the state from subspace 2 is such that all electrons and holes are bound in 
ground state excitons except from one pair which is bound in an excited exciton state. From 
the orthogonality of single exciton wave functions ( |3.4|) with different a it follows that there 
are only two kinds of non vanishing terms. The first are the terms that describe Coulomb 
interaction between different excitons, 

|^K,a(^ei, r hl )\ 2 \ip Kja/ (r e2 , r h2 )\ 2 df el df hl dr e2 df h2 

x [Uee(r e i ~ r e2 ) + u hh (f hl - r h2 ) + u eh (r el - r h2 ) + u eh (f e2 - r hl )] . (3.13) 

This integral converges at large distance, R, between the excitons (dipole-dipole interaction 
falls off as 1/-R 3 ) and it is of the order of (a 2 /S)e}>. The second kind are the terms where the 
identity of the exciton is violated. We gave an example for those terms when we discussed 
the nonorthogonality and saw that they contain the small parameter a 2 /S. Other non- 
vanishing terms, e.g. for states with smaller overlap between subspace 1 and subspace 2 
contain higher powers of a 2 /S. 

From the above arguments we see that the matrix l~i\ exex ) contains the small parameter 
a 2 / S coming from A/12, A/21, ^12 and 7^21- For this reason, in the leading order, all other 
contributions to "j-[( exex ) containing this small parameter can be neglected. In particular, 
in the summation over intermediate states in Eq. (|3.12|) we consider only diagonal terms 



of (7^22 ~~ M 22 E) 1 , and in these terms neglect the interaction between excitons. In other 
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words, TL22 — N22E can be replaced with the diagonal matrix — E, where E^ is the sum 
of the energies of N excitons with at least one of them excited or ionized. 

From this expression the importance of inequality Q3.2| ) is immediately seen. If this 
inequality is not satisfied then the energy of the ground state excitons E can be close to E^ 
due to the high ground state exciton kinetic energy. In such a case some matrix elements of 
^(exex) con £ a j n a sniall denominator and become anomalously large. When condition ( |3.2|) 
is met the exciton kinetic energy can be neglected and E ~ —Neb- As a result, matrix 
elements of "}i^ exex ) are 

n /(exex) 

H -' ~\ E~+m h • (3 ' 14) 

We would like to add one more comment concerning the neglect of the exciton-exciton 
interaction in the intermediate states. The radius of highly excited exciton states can be of 
the order of the distance between different excitons and the interaction between excitons in 
this case is of the order of the electron-hole interaction within the same exciton. However, 
such highly excited states are very close to states of free electrons and holes where Coulomb 
interaction can be neglected compared to the kinetic energy. 

Eq. ( FTTH ) with H {exex) determined by Eq. ((TTJ) is now linear in the energy E, and dif- 



ferent from Schrodinger equation just by the non-diagonal matrix Afn. The transformation 
that reduces Eq. ( p.ll|) to the regular Schrodinger equation is equivalent to the introduction 



of an orthogonal and normalized basis, 

$u = V vv >$ v >, (3.15) 

where the matrix V is not unitary. (Note that this matrix is defined only in subspace 1). 
The corresponding transformation of the M matrix is 

VWV = / , (3.16) 

where I is the unit matrix. Since the difference between M and the unit matrix is small, we 
can write 
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N=I + A, (3.17) 

where A^ = A and A oc a 2 /S. That is, V can be found from Eq. (|3.16| ) with the help 
of the perturbation theory, assuming that V = I + Vi, where Vi ~ a 2 /S. It is necessary 
to note that the transformation ( |3.16| ) is not unique because an orthogonal and normalized 



basis can be chosen by many ways which differ by unitary transformations. So we chose a 
simplest solution to Eq. (|3.16| ) that produces minimal modification of the basis constructed 



of single exciton wave functions, namely, V] = Vi. This immediately gives 

V=l--A. (3.18) 
2 v 1 



After the transformation to the orthogonal basis, Eq. ( |3.11| ) takes the following form 



V ] {Hu - H {exex) )V$>i = E^i . (3.19) 

The matrix Hn can be represented as 

Hu = H +H x , (3.20) 

where Ho describes free excitons and Hi describes their interaction, Hi oc a 2 / S. We neglect 
terms of the second and high order in a 2 /S, so Eq. ( |3.19| ) is reduced to 

H ex ^i = E^i , (3.21a) 

H ex = n n - \ (AHo + HoA) + H iexex) . (3.21b) 



Eq. ( |3.21a|) is the Schrodinger equation for the gas of iV excitons and Eq. ( |3.21b|) is the 



expression for the exciton Hamiltonian that includes exciton-exciton interaction. 

The necessity of the transition to a new basis, Eq. ( |3.15|) , means that exciton-exciton 



interaction changes single exciton wave functions. Such a change appears also in other 
approaches p^j21|j22]] . 



In the leading order in na 2 only matrix elements describing two exciton interaction should 
be kept in Eq. ( |3.21| ). To calculate them it is enough to consider the Hamiltonian of only 



two excitons. This is the subject of the next subsection. 

16 



B. Two-exciton Hamiltonian 



The Hamiltonian of two electrons and two holes has the following form 
fi 2 V 2 , fi 2 Vi ft 2 V?, ft' 2 V?, ,,,,,, 

where 

u (r el ,r hl ,r e2 ,r h2 ) = u ee (\f el - r e2 |) + u eh (\f el - r hl \) 

+u eh (\r e i - r h2 \) + u eh (\r e2 - r h \\) + u eh (\f e2 - r h2 \) + u hh (\f hl - f h2 \) . (3.23) 

For the calculation of matrix elements of Tin and Afn it is necessary to know only the 
wave function (3"!3) of two ground state excitons. Explicitly, this function is [[J 

®i? u s r ,K 2 ,s 2 (rel<Tel,rhl<7hl,re2Ve2,r h 2Vh2^ = ^ 

X [® Kus^el, (Tel] f hl , a hX ) ^ ^ (f e2 , O e2 \ f h2 , (T h2 ) 

-^K ltSl if e2 , a e2 ; f hl , cT M )^ /?2 S2 (f el , a el ; f h2 , a h2 ) 
- *^i,«i(*Vi> r h2 , (T h2 )^ g 2 S2 (f e2 , a e2 ; f hl , a hl ) 

+^K uSl (re2,(Te2;r h2 ,a h2 )^g 2S2 (r el ,a e i;r hl ,a hl )\ . (3.24) 

Hereafter the quantum number a, characterizing an internal exciton state, will be omitted 
in the case of the ground state. The unity matrix elements in the space of two ground state 
excitons are 

^K 1Sl ,K 2 s 2 ;K 3 s 3 ,K4S 4 = ^ S ^ S 3^ S 2S4.^K 1: K 3 ^K 2: K 4 + ^s 1 s^s 2 S3^K u K 4 ^K 3 J<4 • (3.25) 

The two products of 5-symbols appear because the states are symmetric with respect to the 
transposition of single exciton quantum numbers. 

The matrix An in the subspace of two excitons states is J^fg lS k 2 s 2 -k 3 s 3 k 4 s 4 = 



{^Ki s V k 2 s 2 \^k 3 s 3 -k 4 si)- According to the two exciton wave functions definition, ( 3.24| ), 
Eq. (^T7|) has now the form ^k iSi ,k 2 s 2 ;K 3 s 3 ,k 4 s 4 = ^^k^k^k^ + ^k iSi ,k 2 s 2 ;K 3 s 3 ,k^ 
where 
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J ^K 1 s l ,K 2 s 2 ;K 3 s 3 ,K 4 s 4 g ^S 1 +K 2 ,K 3 +K 4 

X {^Sl e S 4e ^S 2e S 3e ^S lh S 3h ^S 2h S 4h ^«leS3e ^«2e S4e ^lh S 4h ^ S 2h S 3 h ) ^ ' (3.26) 

s je and SjTj are electron and hole spins respectively of the exciton with the spin Sj, and 
where 

<f> q = [ e"^0(r) df , (3.28) 



is the Fourier transform of the wave function. The reduction of the overlap integral to such 
a simple form is possible due to the small parameter, Eq. ( |3.2| ). Details of the calculation 
are given in Appendix |A 1| . 

For the calculation of the matrix elements of Tin it is convenient to separate them into 
two terms, 

(Tin) - - - - = 4- CS 29") 

V L H/ K lSl ,K 2 s 2 ;K 3 s 3 ,K 4 s 4 K lSl ,K 2 s 2 ;K 3 s 3 ,K 4 s 4 K 1 s 1 ,K 2 s 2 ;K 3 s 3 ,K 4 s 4 ' \ ' J 

The first part, H^ d \ contains matrix elements between initial and final states where excitons 
consist of the same particles. This part consist of 8 integrals which have similar form. The 
grouping of terms which are different only by the notation of integration and summation 
variables leads to 

H (d) 



K 1 s 1 ,K 2 s 2 ;K 3 s 3 ,K 4 s 4 

0"elO"e20"/il07i 2 ' 



J df el df hl df e2 df h2 (3.30) 

s2&hl a 'h2 

X^ isi (r e icr e i, r W CT W )*£ 2S2 (r e2 cr e 2, r h2 a h2 )H 2e h^ K SS3 (r el a el , r hl <7 hl )^ £ 4S4 (r e2 a e2 , r h2 a h2 ) 
+ J df el df hl df e2 df h2 . 

o- e lo- e2 o- hl a h2 

X^*g isi (f el a e i,f hl a h i)^f*g 2S2 (f e2 a e2 , r h2 a h2 )H 2eh ^ ^ 4S4 (f el a el , r hl a hl )ty g 3S3 (f e2 a e2 , r h2 a h2 ). 

The second part, H^ x \ contains matrix elements between initial and final states where 
excitons consist of different particles. Again we have 8 integrals that can be grouped to the 
following form 



- - 

E"lSi,_ftr 2 S2;^3S3,^'4S4 



J df el dr hl df e2 df h2 (3.31) 

&elO'e20'hl< 7 h2 

x ^*k 1Si (reiVei, r hl a hl ) ^*g 2S2 (r e2 o e2 , r h2 a h2 )H 2eh ^ g jSs (r e2 a e2 , r h ia hl )^g 4S4 (f el a e i,f h2 a h2 ) 



X 



I df el dr hl df e2 df h2 . 

$g isi (r el a el , r hl a hl )^/g 2S2 (f e2 a e2 , r h2 a h2 )H 2eh ^ g 4S4 (r el a el , f h2 a h2 )^ g 3S3 (f e2 a e2 , r hx a hl ). 



There is a confusion in the literature about the term "exchange" in the description 
of interaction (see, however, Ref. |J, where clear definitions are given). For elementary 
particles, the first term in Eq. ( |3.30|) is usually called " Hartree interaction" and the second 
is called "exchange interaction" while H^ x ' does not exist. In the case of excitons is 
called "direct interaction" and H^ x ' is called "exchange". We stick to this last terminology. 

In the direct part, ( |3.30| ), there are terms of H, 2e h which do not mix between the two 
excitons, i.e., the kinetic energy of all the particles and the potential energy between particles 
belonging to the same exciton. These terms will result with the unity operator (|3.25|) 
multiplied by the single exciton energy. The other parts of 7i 2e h are the Coulomb interaction 
between the two excitons. The details of the calculation are given in Appendix [A2| and the 
result is 



H 



{d) (3.32) 



where Ek = — e& + h 2 K 2 /2M is the single exciton energy, and 



U 



id) 



Kisi,K 2 S2,K 3 s i ,K4S4 



6 SlS3 6 S2Si U ! ^ d \\K 1 - K 3 \) + 6 SlSi 6 S2Sa U (d \\K 1 - K 4 \) . (3.33) 



Here 



U^(q)=u ee (q) [</> 



qm h /M 



+ u hh (q) (j) 



qm e /M 



+ 2u eh {q) (f) 







;< 






qm h /M 




qm e /M 



where Uij(q) is the Fourier component of the Coulomb potential, 



Uij(q) = j Uij(r)e %qr dr 



(3.34) 



(3.35) 
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and [(f) ] is the Fourier transform of the wave function squared, according to the general 
definition, 

[(jxj> a ) q = J e-^4>{r)(j> a {r)dr . (3.36) 

In the exchange part of the Hamiltonian matrix element, ( |3.31| ), all the terms mix be- 
tween the excitons since the initial and final state consist of different particles and therefore 
corresponds to different pairs of excitons. The single exciton energy is now multiplied 
by the factor A which reflects the small overlap between the initial and final states. The 
result of the calculation is 

tfM _ U(X) - 4 ^ s 

X {fis lt ,S4 se fis2 f ,S 3e fisi h S 3h fis2hS4h ^Sl e S 3e ^S 2e S4e ^S lh S 4h ^S 2h S 3h ) ) (3.37) 

where A was defined in Eq. (|3.27 ) and 

^ ' r r ... ,/,7 I 2 

dr . (3.38) 



[,{ "^/* J ''(|-/ Wr)+ "" wl 



l2 iqr d Q 

q (2tt) 2 



Details of the calculation are given in Appendix |A3 . 



Apparently, the first term in H# T -> r -t r > , Eq. (B.32j), describes free excitons. 
According to the notation of Eq. (|3.20| ) this is a matrix element of H.q. In other words, 

H-0K 1 s 1 K2S2,K 3 s 3 ,K 4 s 4 , = {^s 1 s 3 &S2Sa&K 1 ,K 3 ^K 2 ,Ka ^ s ^ s ^ s ^ s ^Kx,K^R 2 ,R 3 ) (Er 3 + ^K A ) ' (3-39) 

The second term in r -> - r> , as well as H^§) r> r , r> , describe exciton- 

K 1 s 1 ,K2S2,K 3 s 3 ,K 4 S4' K 1 s 1 ,K2S2,K 3 s 3 ,K 4 S4' 

exciton interaction and, according to the same notation, are parts of the matrix element 
of Hi. 

Matrix elements Ti.^, seem to be of higher order in the small parameter, a 2 / S, since it 
contains the product of two off-diagonal matrix elements, (H^+A^A^e^^T^i+A^iAfe;,)^/, 
each of which is proportional to a 2 /S. However, the additional factor of a 2 /S is cancelled 
by the summation with respect to the momentum of intermediate states in Eq. ( |3.14 ). 



Therefore the order of magnitude of TC UU , is the same as of other matrix elements describing 
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pair exciton-exciton interaction, namely, e&a 2 / S. Similar contribution to the exciton energy 
appears also in other approaches where it is usually referred to as screening correction. 
P J2T|J30|J3T| The calculation of Ti.^, ex \ is much more complicated because they contain wave 



functions of excited exciton states. The details of the calculation of Ti vv i are given in 
Appendix |A4|. The result has the form 



H 



(exex) 



I A 

g °K 1 +K2,K 3 +K4 



where 



(^S e lS e 3^S hl S h 4^S e 2S a 4Ss h 2S h3 + ^S el S e 4^S hl S h3 Ss e2 Se3^ 



H2 s h4. 



v 2 



Vl 2 / E ai + E a2 + h 2 q yM + 2e b (2tt) 2 ' 



(3.40) 

(3.41a) 
(3.41b) 



,„.„,., E ai + E a2 + h 2 q 2 /M + 2e b (2vr) 2 ' 

Here in the sums with respect to the internal exciton quantum number a at least one of cxi 
and a 2 corresponds to an excited state. The remaining matrix elements are 

D< &*M = U M [#ai]-m||g7M [#Jr»||<f/M + U hh{q) [WaA m ^/ M [#<* 2 ] - m ^/M ( 3 - 42a ) 
+U eh (q) [#ai]_ m , |S/J if [00a 2 ]- me?7 M + U eh(q) [#ai] m „#Af [#a 2 ] meg7 M > 



^Sjtf) = / [0a 1 ,-k2+qrn l] /M0a2,-k 1 -q m]l /M U ee(\kl ~ &2 + q\) 



(3.42b) 



dk\dk 2 



+ 



K,-ki-qrn e /M0a2,-k2+<rrn e /M U hh(\kl ~ k 2 + q])^^-^ 
U eh (\kl - k 2 \) (0fc 2 </>|fc 2 _ 9l + hx^-q]) 



ai,fci— qrrie/MT a.2,k2— qm^ /M 



dk\dk 2 



According to Eq. ( |3.21b|) , to complete the calculation of the interaction Hamiltonian it 



is necessary to calculate AHq + HqA. Both A 



V\V2,VZV4 



and TC, 



§V\V2,VZV4 



are symmetric with 



respect to transpositions of the first and the second pair of subscripts, i.e., 



A = A = A 

• / ^-V\V2,V?,V4 J ^-V2V\,V- 3 ,V4 • n -V\V2,V4V% 1 



7~Lc 



Qv\V2,vi,V4 T~~tov2V\,V3V4 riQv\ V2,V4V% 



(3.43a) 
(3.43b) 
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They are defined in the symmetric space, i.e., (1/11/2) and {v^vi) describe the same two- 
exciton state. Because of this, the summation with respect to intermediate states in 
A Ul i, 2tU3 u 4 ,'Hou 3 v 4 „u 5 u 6 has to be carried out over all different pairs of z/3 z/4, but not indepen- 
dently with respect to z/3 and z/4. Actually, due to continuous spectrum of K the contribution 
of diagonal states, where z/ 3 = z/ 4 , can be neglected and then the independent summation 
respect to z/3 and 1/4 is equivalent to double counting of each pair Z/3Z/4. As a result it is 
possible to sum independently with respect to z/3 and z/4 and then to divide the result by 2. 
This leads to 

(x) A 
{AHo + W ^)| iaiAsaA , 3Aa4 = ^ $k 1+ k 2 ,k 3+ k 4 (3-44) 

X \^S le S 4e ^S 2e S 3e ^S lh S 3h ^S 2h S 4h + S Sle S 3e 8 S 2e «4e ^ S lh S 4h ^ S 2h S 3h ) 



4Ae b 



X (^Sl,;S4 e ^S2 e S3 e ^S lh S3 h 5s2hS4h "I" ^Sl e S3 e < 5s2eS4 e ^lhS4 h ^S2hS 3 , l ) • 



Here we used again the assumption of a small exciton kinetic energy, ( p.2| ). 

Finally, we can present matrix elements of the two-exciton Hamiltonian in the form 

H 2 ex = ftwis u it2S2,K 3 s 3 ,K A S4, + g ( U K lSl ,K 2 s 2 ,K 3 s 3 ,K 4 s 4 + U ^ lS ^K 2 s 2 ,K 4 s 4 ,K 3 s 3 ) ' (3.45) 



where matrix elements of 7io are defined by Eq. (|3.39|) and matrix elements of U ex ex are 



u z7:k 2S2 ,k 3S3As4 = - k 3 \) + v x \ 

+ 6 K 1 +K 2 ,K 3 +K i 6 s le s 4 Js 2e s 3 Js lh s 3 Js 2h s 4h (U {x) - 2Ae b - V 2 ) . (3.46) 

In this form we still keep a trace of non-elementary nature of excitons in the spin 5-symbols. 
The corresponding expression in exciton spins is quite cumbersome. We overcome this 
disadvantage in the many-exciton Hamiltonian in the next subsection. 
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C. Hamiltonian of exciton gas with pair interaction 



The matrix form of the Hamiltonian of the exciton gas ( |3.21b| ) is practically inconvenient. 
It is desirable to reduce it to the second quantized form. There is a standard way to obtain 
the second quantized form of the Hamiltonian of a gas of elementary bosons. To make use of 
this way we show that the Hamiltonian of two excitons (|3.45|) is equivalent to a Hamiltonian 
of two elementary bosons. Then, keeping only pair exciton interaction, we can immediately 
write down the Hamiltonian of iV excitons. 

The Hamiltonian of two elementary bosons with masses M and the interaction U can be 
written as 

h 2 h 2 

H ex = -2e b - 4t V? - 4t + U , (3.47) 
2M 2M 

where e& is the energy necessary to create a boson. 

As a basis for wave functions we chose plane waves so that wave function of two bosons 
can be written as 

1 



Kisi,K2Sa ^y2S 



g Sl (a 1 )g S2 (a 2 )e i ^ +K ^ + g^Mg^Me*********) , (3.48a) 



when at least one of the inequalities S\ ^ s 2 , K\ 7^ K 2 is satisfied and 



= s > ( 3 - 48b ) 



when both spins and wave vectors of the bosons are equal. 

Matrix elements of the Hamiltonian, ( 13.47 ), between functions, ( |3.48a|) , are identical to 



the matrix elements, ( 3.45Q , if we identify the matrix elements of U between wave functions, 



(|3 48a|) U e -?~ ex ^ ^ - -\-U e *~ ex ^ - - with jj ex ~ ex _ _ \jjex-ex 

VtL ij/' K 1 si,K2S2,KsS3,K4S4 KlSi,K2S2,K4S4,K3S3 KiSi,K2S2,KsS3,K4S4 KiSi,K2S2,K4S4,KsS3 ' 

Eqs. (|3.45| )-( p.46| ). This last identification is the definition of U. The nontrivial spin 



structure of these matrix elements means that the operator U is spin-dependent. [One 
should remember that Sj = Sj e + Sjh, = dj e + cr^,, and g s (a) = g S( ,(a e )g Sh (ah)] Due to 
the continuity of K, the contribution of matrix elements with K\ = K 2 and with K 3 = K 4 
to all effects is negligible and we don't consider them. 
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Now we have two-exciton Hamiltonian with known matrix elements of the pair interaction 
operator. The many excitons Hamiltonian with pair interaction in the second quantized form 
is usually written with the help of a non-symmetrized matrix element, f3l 



X 



J dR x dR 2 e ^<*^ + *^>tf e ^ M e i ^+*^ 9s a ((Ti)gs.M ■ (3.49) 



Due to commutation relations of Bose operators it is possible to write down this Hamiltonian 
also with the symmetrized matrix element, (l/2){U e ^~ ex J> - - + U^' e % r> - ). 

K\S\,K2S2,K 3 S 3 ,K4S4 KlSi,K2S2,K4S4,,KsS3 / 

That is, an additional factor of 1/2 appears in the interaction term of the many exciton 
Hamiltonian 

Hex = e kc ] K)S c K)S + (3.50) 

K,s 

+7o E {[{U^m -K 3 \) + V 1 ) 8 SUS3 S S2 , S4 + (uWQlb - K A \) + V,) 5 Sl , s J S2 , S3 

3 1 ,S2' S 3< 3 <i 

+ Kc (^S e l,S e 4^S e 2,S e 3^S hl ,S h3 ^S h2 ,S h4 + ^Sel,S e 3^S e 2,S e 4,^S h l,S M ^S h2 ,S h3 JJ 
XQ K 1 +K2,K 3 +K4 C K2,S2 C K 1 ,S 1 C K 3 ,S 3 CK4,S4- 

where cg s and are exciton annihilation and creation operators respectively, and 

V x = U (x) - 2Ae b - V 2 . (3.51) 

Now we make use of the one to one correspondence between electron and hole spins and 
the exciton spin, in order to express the ^-symbols containing separately electron and hole 
spins, in 5-symbols containing exciton spin only. Details of this calculation are given in 
Appendix [B], and the result is 

H ex = E K C ] KyS C KtS + — [ Uidi {Q) + ^l] c t a ,-a C k,*i C «i-9.»i C ^2+9,^ 



25 ^ 

K,S K 1 ,K 2 ,q 

S 1. 3 2 



v x x , 

Ki,K 2 ,q 



K\s\ C K 2 -s 1 C Ki -qs 2 C K 2 +q-s 2 

L S1S2 



~^Y cJ< K ls ( ^K 2 - s C K 1 -qs C K2+q-s + 2 E C i? lSl C i? 2S2 C Xi-gsi C i? 2 +gs 2 

S S1S2 
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(3.52) 



Hamiltonian, ( |3.52|) , is the main result of this section. In the following sections we use 
this Hamiltonian to study the density dependent luminescence line shift and exciton-exciton 
relaxation. To do this we calculate numerically the matrix elements of the Hamiltonian for a 
number of structures. The calculation of V\ and V2 is extremely difficult because it involves 
all exciton excited states. Fortunately, for the most of the structures that we consider V\ 
and V2 are numerically small compared to other matrix elements. The reason behind this 
is that the wave functions of excited states oscillate in the region where the wave function 
of the ground state is smooth. Our estimates show that the value of 7i^, ex ^ is smaller then 
10% of Ae b or U^ x \ So in the following calculation V\ and V2 are neglected. 



IV. MEAN FIELD APPROXIMATION. LUMINESCENCE LINE SHIFT 

In the mean field approximation the exciton scattering is neglected and only the shift 
of the exciton energy due to interaction is taken into account. This means that only terms 
diagonal with respect to occupation number ng a = ct cg g are kept in the Hamiltonian. 
Then the energy of the system of N excitons is 

E N = Y: E K n Rs + ^ 2 + ^|E (N. - N_ s f 

+ £ U^m-K^n^n^ (4.1) 

K u K 2 ,s 

where N s = J2k u k s * s num ber of excitons with spin s, N = J2 S N s , and 

V b = U^(0) + V x . (4.2) 

The recombination energy eg of an exciton with momentum K and spin s equals the 
change of the energy of the exciton system when the occupation number n^ s decreases by 
one, eg s = E^{ng s ) — Ejsr{ng s — 1) ~ dE^ / dng s . Neglecting the photon wave vector we 
can put K — and then Eq. Q4.1|) leads to 

e s = -e b + V b n + V x (n s - n_ s ) + J ^U^(K)ng s . (4.3) 
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where n s = N s /S is the density of excitons with spin s and n = N/S is the total exciton 
density. The energy splitting between optically active excitons with spins +1 and -1 is 

A-E^-i = ei - e_i = 2T4(m - n_i) + | 0^U (d) (K)( nii l - ng ^) . (4.4) 

Eqs. ( |4.3D and ( |4.4j ) are complicated to use due to the integral term. To calculate it, 
the exciton concentration is not enough, it is necessary to know the energy distribution of 
excitons. These expressions can be simplified if the typical exciton wave vector is much 
smaller than the typical wave vector of U (d) {K). The last one is characterized by the wells 
width, L and the width of the barrier, w (see Appendix The exciton wave vector is 
small in the case of resonant pumping at low temperature T. Then it is of the order of the 



photon wave vector K p h or of the thermal wave vector ~ y/2MT/h. So, if both K p h and 
Kt are much smaller than 1/L and 1/w then U^ d >[K) in Eq. ( f4.1|) , (|4.3| ) and ( f4.4|) can be 
replaced with U^(0). As a result we have 

E N = Y, E K n Rs + ^N> + ^ (N s - AL S ) 2 + ^ £ iV s 2 , (4.5) 

and 

e s = -e b + V b n + V x (n s - n_ s ) + U {d) (0)n s , (4.6a) 

A£ w = Vesim - n-i) , (4.6b) 

where V es = 2V X + U^(0). An expression similar to Eq. (|4.6a|) was suggested phenomeno- 
logically by Amand et al. This expression differs from Eq. ( |4.6a|) by a relation between 
the coefficients in two last terms. 

From Eq. ( |4.5| ) we see that for a constant exciton number, N, the ground state of the 
exciton gas can be paramagnetic or ferromagnetic, depending on the sign of V es . When V es is 
positive the minimal total energy of the system is reached when the number of excitons with 
opposite spins is equal, iV s = iV_ s . For negative V es , the system reaches its minimal energy 
when the difference, N s — AL S , is maximal, which corresponds to ferromagnetic phase. The 
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possibility of these two phases for condensed excitons has been pointed out by Fernandez- 
Rossier and Tejedor. [[33|] From Eq. ( |4.6b| ) we see that the same parameter, V es , characterizes 
the spin energy splitting. 

To make a quantitative comparison to experiments we evaluate the values of the Hamil- 
tonian matrix elements for coupled quantum wells where electrons and holes are confined in 
different wells. Such a separation is usually reached by an external electric field applied in 
the growth direction. In the region of well widths that we consider the energy separation 
between the ground state and the first excited state in the wells is much larger than typical 
external potential drop across one well. This means that the electron and hole ground state 
wave functions can be taken as 

r~2~ nz e 

(e(z e ) = J — sin -j^- , -L e < z e < , (4.7a) 

. , . r~2~ . ir(z h — w) . . 

C,h(z h ) = \ — sin , w<z h <w + L h , (4.7b) 

V L h J^h 

where L e and and the widths of the electron and hole quantum wells respectively and w 
is the width of the barrier. The expressions for the interaction energy Uij are cumbersome 
and their Fourier transforms are presented in Appendix 0. 

For the calculation of matrix elements we use the variational single-exciton wave function 
that gives a very good approximation |23], 



l2ub{b + r ) I 2b 



where b and r are variation parameters that are found by minimizing of the binding energy. 
The Fourier transform <p q necessary for the calculation of (q) and V x is given in Appendix 
0. 

For Z/W(0) it is possible to obtain a simple analytic expression without making use of 
the exciton wave function (Appendix |C|) , 

U(d) (o) = [ w + 0.397(L e + L h )}. (4.9) 
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We see that the direct interaction, U^(0), grows with the separation between the wells, d = 
w+(L e +Lf l )/2. This behavior is easy to understand. From the point of view of electrostatics, 
excitons resemble parallel dipoles of the size d, and U^ d > (0) is a the dipole - dipole interaction. 
This interaction is a repulsion growing with the size of the dipoles. The coefficient of the 
second term in the square brackets is close to 0.5, for which Z7^(0) corresponds to the plate 
capacitor approximation P, 51|,P5] . Typical value for £/w (o) for coupled quantum wells where 
the separation between the wells is d — 100A is about 1.5 x 10~ 10 meV- cm 2 . 

We present numerically calculated V x , V& and V es for a symmetric ALjGai^As / GaAs 
/ Al x Gai_ x As / GaAs / Al^Gai-^As coupled quantum wells structure with the well widths 
L e = Lh = L. In the numerical calculations we use the electron effective mass m e = 0.067mo, 
the hole effective masses m± = 0.45mo and mil = 0.126mo (here mo is the free electron mass), 
and the dielectric constant k = 12.5. In the previous paper [^J we have shown that the 
exciton wave function depends mainly on the distance between the centers of the wells and is 
not very sensitive to details of their geometry. We can expect the same from the parameters 
V x , Vb and V es . For this reason in Figs. [1] - f| we present the dependence of V x , Vf, and V es on 
this distance, d = w + L. To show the sensitivity of these parameters to the barrier width 
and the well width separately we give on Figs. |] and |3] two curves, one for a given w (solid 
line) and the other for a given L (dotted line). 

The parameter V x can be positive or negative depending on the separation between the 
wells d. To demonstrate this we presented its dependence on d for different ranges in Figs. 
|TJ and 0. The reason for the change of the sign is that V x contains two contributions of 
different sign, the negative quantity and the positive quantity — 2Aeb, Eq. ( [3.51|) (we 
neglect V2). At small distances between the wells the second term dominates. When the 
separation between the wells increases U^ x ' is just weakly affected while the binding energy 
6b decreases, which leads to the change of the sign. 

In Fig. m we present the parameter Vb which characterizing the overall shift of the exciton 
luminescence line. It is the combination of the direct interaction and the exchange term, V x , 
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Eq. (|4.2| ). The direct interaction dominates and V& is always positive, which leads to the 
blue shift of the line. 

To demonstrate the possibility of paramagnetic and ferromagnetic phases of the exciton 
gas we present the dependence of V es on d in Fig. |. We see that V es is a decreasing 
function of the separation. It becomes negative at large separation, which corresponds to 
ferromagnetic phase. 

In a single infinite quantum well the electron and hole wave functions describing the 
confinement in the well are equal, and according to Eq. (|3.34j ) U^(q) is identically zero. 
HU For an estimate of the exchange matrix element V x in a single-well we can use two- 
dimensional model for which the exciton wave function is known (see, e.g., Ref. 0) and 
which is typically used for such estimates. [|],[7|-|3 In this case the wave function is a simple 
exponent, the absolute value of the binding energy is e& = 2/ie 4 /fr 2 K 2 and 

e b A = -— — , (4.10a) 

V. = ^ (l - » 3.03 f l . (4,0b, 

/i \ 4096 J fj, 

The corresponding coefficient in the expression for the splitting, ( 4.6b| ) appears to be twice 
larger than in Ref. || . This difference comes from a different numerical factor in the Hamil- 
tonian. 



V. EXCITON-EXCITON RELAXATION TIME. 

In this section we study exciton gas relaxation due to exciton-exciton collisions. In 
these collisions, the excitons change their momenta and can change their spins, however, 
the sum of exciton spins is conserved. For further calculation it is convenient to introduce 
the notation M^^(q) for the scattering matrix element that describes the scattering from 
the state with spins s± and s 2 to the state with spins s 3 and s 4 (si + s 2 = s 3 + s 4 ) with the 
transferred momentum Hq. In the Born approximation, according to Hamiltonian ( |3.52| ) the 
spin dependence of the matrix element is reduced to the separation of three cases. If the 
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sum of exciton spins is nonzero then the collision matrix element 



M s s i:i{q) = U {d \q) + V x 



s 2 + st ^ 



(5.1a) 



If the total spin is zero and spins in the final state are different from spins in the initial state 
then 



S2^ S X . 



(5.1b) 



If the total spin is zero and initial and final spins are the same then 

M-(g) = ^(<?)-y • 



(5.1c) 



For 2D scattering, however, the Born approximation at low energies is not satisfactory 
and the matrix element can be strongly renormalized. P^-B3 For a spin independent inter- 
action between particles, U(r), the renormalization is reduced to the division of the Born 
matrix element by 



M 



1 1 /_ . Kr 

h - C + In 

2i 7T V 2 



U(r)df 



(5.2) 



where h 2 K 2 /M is the kinetic energy in the center of mass reference frame and C is the Euler 
constant. We don't calculate the renormalization for a spin dependent potential since we 
use Eq. ( |5.2j ) only for estimates. 

The spin relaxation in the exciton gas is usually described by simplified kinetic equations 
that ignore exciton momentum distribution. |T7| , |39| . The Bolzmann equation that describes 
both spin and momentum relaxation has the form 

f)TI -* 27T C 2 r 



dt 



-( n K + qs + !) + 1 ) n Ks U K 1 s 1 

d 2 K x d 2 q 



x6 [E K + E Kl - E^-E^j (2?r)2 (2?r)2 



2vr 



+ y E / \M s s r s ^(q) K s + l)(n Si _ s + 1 
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X 



d 2 K x d 2 q 



S (e k + E Kl - E R+? - E Rx _ 4 



(5.3) 



(2tt) 2 (2tt) 2 ' 

Numerical solution of Eq. ( |5.3| ) is much more difficult than its simplified versions that 
ignore momentum distribution. In general, the relaxation according to Eq. ( |5.3j ) cannot be 
exactly described by a relaxation time. To characterize the relaxation rate it is possible, 
nevertheless, to introduce an inverse relaxation time as the coefficient for —ng s in the 
collision operator 

1 2vr 



h 



E / M "l (?) [( n K+qs + n K 1 ^gs 1 + 1 ) Tl K 1 s 1 ~ n K+gs n K 1 -gs 1 



27T 



x5(E K + E Kl -E^-E 
2 



d 2 K x d 2 q 



K+g K.-qJ ( 2?r )2 ( 2 7t) 2 



y n K+q Sl + n K 1 -q-s 1 + ^) n Ki-s 71 K+q Sl n K 1 -q-s 1 



xS[E K + E Kl -E 



K+q 



E K^q) 



d 2 K x d 2 q 



(5.4) 



(2tt) 2 (2tt) 2 

This expression contains both linear and quadratic terms in exciton occupation numbers. 
Respectively, an order of magnitude estimate of the relaxation time contains linear and 
quadratic terms in the exciton concentration n, 

1 n 2 + nK 2 



MtlliK) 



(5.5) 



t Ks h(h 2 K 2 /fi) 

The relaxation time tk s characterizes energy and momentum relaxation in the exciton 
gas. The same time characterizes also a partial spin relaxation. The spin relaxation due to 
collisions cannot be complete because of the total spin conservation in collisions. Formally 
it is described by the identity 





Of 



(n s - n_ s ) = 



(5.6) 



that follows from Eq. ( |5.3| ). Here n s is the concentration of excitons with spin s. That is, the 
only spin relaxation due to collisions is the relaxation between dark excitons, (excitons with 
spin ±2), and bright excitons, (excitons with spin ±1) (see also Ref. ||). It immediately 
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follows from here that spin relaxation in the exciton gas is characterized by few relaxation 
times which correspond to different relaxation mechanisms. If the exciton concentration is 
not very small then the fastest relaxation is the dark - bright exciton relaxation characterized 
by Tks- Complete relaxation can take place due to processes that involve D'yakonov - Perel 
mechanism of electron spin relaxation, light - heavy hole mixing, or electron - hole exchange. 
p5| All these mechanisms contain their respective small coupling constants and are activated 



by scattering. So, if the main exciton scattering mechanism is exciton-exciton scattering then 
the complete spin relaxation has to be much slower than dark - bright exciton relaxation. The 
relation between the relaxation times can be different for a very small exciton concentration 
when other scattering mechanisms, e.g., phonon, impurity or surface roughness scattering 
are important. 

VI. DISCUSSION 

In this section we compare our results with a few experiments. We consider the shift of 
the exciton luminescence line, the energy splitting between exciton with different spins, the 
polarization of the exciton gas, and the exciton-exciton scattering time. 

As we already mentioned in the introduction, the density dependence of the exciton 
luminescence line shift, the time dependence of line spin splitting and the luminescence 
depolarization fl0|-13, r6l,|l7| proved that all these phenomena come from exciton-exciton 



interaction. A theoretical study of these phenomena has been done by Fernandez- Rossier et 



al. [§,^11, Ciuti et al. ||, Amand et al. ||17|| , and Maialle et al. ||39|| . Here we consider only 
some features of the experiments which have not found a clear explanation so far. 

The density dependent blue shift of the exciton line in GaAs/Al z Gai_ x As symmetric 
coupled quantum wells, where electrons and holes are spatially separated by external gate 
voltage, has been recently measured by Butov et al. |J and Negoita et al. ||. Butov 
et al. detected blue shift of 1.6 meV at zero magnetic field for wells width L = 80 A 
and barrier width of w = 40 A. They used the plate capacitor expression for the direct 
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interaction, neglecting the exchange, to calculate the exciton density of n = 9 x 10 9 cm -2 . 
To compare our results with their measurements we calculate the concentration according 
to Eq. ( |4.3| ) and check the importance of the corrections. For the described geometry we 
obtain V b = 8.7 x 10" 11 meV- cm . The temperature of this experiment, 50 mK, is so low 
that we can use Eq. ( |Q| ) instead of Eq. (fO|), where 17^(0) = 1.5 x 10~ 10 meV- cm -2 . 
Assuming equal concentration of excitons with different spins, n s = n_ s = n/4, we have 
5e s = Vbn + (0)n/4 and for 5e s = 1.6 meV we obtain n = 1.3 x 10 10 cm -2 , which is 
close to the concentration obtained from the simple plate capacitor expression without the 
exchange correction. 

The comparison with the results of Negoita et al. || is more interesting because there 
the excitation concentration was measured from the excitation intensity, independently of 
the blue shift. The measurements were made in symmetric coupled quantum wells, where 
the wells width was L = 60 A and the barrier width was w = 42 A. The lattice temperature 
was 2 K and the concentration was in the range 10 9 — 10 12 cm" 2 . For low density, linear 
blue shift of 5 x 10~ n meV • cm 2 , was observed. For this geometry we have V& = 7.8 x 10~ n 
meV • cm 2 . In the case, n s = n_ s = n/4, with U^(0) = 1.3 x 10~ 10 meV- cm -2 we get from 
Eq. ( }4.6| ) linear shift of 1.1 x 10 -10 meV • cm 2 . The difference between this value and the 
experimental one can result from our assumption that electrons and holes are completely 
confined in separate wells which increases the direct interaction . Another reason can be 
the presence of free carriers in the experiment, which screen the Coulomb potential and make 
exciton-exciton interaction weaker. The later possibility is supported by the luminescence 
line width which is larger than the exciton binding energy. 



In the same structure, in a weak magnetic field, Snoke et al. |4(| observed a red shift 
of the exciton luminescence line which grew with the gate voltage that separated electrons 
and holes. The most striking result is that the red shift reaches values of 10 or 20 meV 
(depending on the gate voltage) at magnetic field around 1 T. Such a magnetic field is not 
strong enough to induce a significant blue diamagnetic line shift. pi|,|4"2"|] The explanation 



that we suggest is based on a very narrow line of the pumping laser. Even a weak magnetic 
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field can split the exciton lines with different polarizations so much that they go away from 
the resonance with the pumping laser. This leads to a reduction of the exciton density, 
resulting, according to Eq. (|4.6|) , with a red shift of the luminescence line. An increase of 
the external electric field increases the separation between electrons and holes leading to 
larger values of the coefficients Vj, and U (0). As a result, the red shift also increases, as it is 
observed in the experiment. The effect is symmetric to the direction of the magnetic field, 
which is also in agreement with the experiment. The magnitude of the effect depends only 
on the absolute concentration change, An, so the relative change An/n in the luminescence 
intensity, can be small. 

Another phenomenon related to exciton-exciton interaction is a spontaneous energy split- 
ting between excitons with opposite spins. A typical experimental way to produce a polarized 
exciton gas, (i.e., a gas where n_ s ^ n s ) is pumping by polarized light. In Refs. [[H|-[l3|, the 
energy spin splitting was measured in multiple quantum wells where the electrons and holes 
are in the same well. This corresponds to zero separation between the carriers and therefore 
zero U^(0) and positive V es . In all the experiments the spin majority excitons had higher 
energy than the minority and the difference increased with the density, as we would expect 
from positive V es . Another evidence for V es being positive in multiple quantum well systems 
is the depolarization of the initially polarized exciton gas, that was reported by different 
authors H JT0|JT^Jl3|jr^| . According to the results of Sec. |V|, when V es is positive the system 



is paramagnetic, and the minimal energy of the system is reached when the exciton gas is 
depolarized. 



In a double well structure Aichmayr et al. |14| detected an energy spin splitting depen- 
dence on the gate voltage that separated electrons and holes. As the voltage increased, the 
energy splitting decreased from 4 meV to zero. This behavior corresponds to Eq. (4.6b|), 



where with the increase of the electron and the hole separation, the coefficient V es , being 
positive, decreases to zero (see Fig. A more detailed consideration, with the help of Eq. 
(4.6a), can describe a different behaviour of the minority and majority exciton luminescence 



lines which is presented in Fig.2 of Ref. JL4J]. At low gate voltage the lines are split nearly 
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symmetrically (the shift of the majority line is positive while the shift of the minority line 
is negative) with respect to the value to which both of them relax with the time constant 
r s d = 180 ps. With increase of the voltage the shift of the majority line does not change 
while the shift of the minority line decreases, becomes positive and at a very high voltage the 
splitting disappears. First of all it is necessary to note that the luminescence decay time (400 
ps for the low voltage and 1000 ps for the high voltage) is a few times larger than r a d- The 
luminescence decay characterizes the decrease of the exciton concentration and the compar- 
atively small value of r s d means that the reason the splitting relaxation is not decrease of the 
exciton concentration but spin relaxation, probably due to light and heavy hole mixing. So 
the line shift due to exciton polarization under the condition of constant total concentration 
can be calculated according to Ae s = V x (n s — n_ s ) + U^(0)n s . At low gate voltage the 
separation between electrons and holes is small, the direct interaction U^(Q) is negligible 
and the shifts of the majority (+1) and minority (-1) lines are Ae±i = ±V x (n + i — 
This is a symmetric shift in agreement with the experiment. With increase of the voltage 
the separation between electrons and holes grows leading to growth of U^(0) and decrease 
of V x . Given the concentrations, at some intermediate voltage V x = (0)n_i/(n + i — n_x) 
and then Ae + i = f/^(0)(n + i + n_ 1 ) while Ae_i = 0. This corresponds to Fig.2b of Ref. |TJJ. 
At high voltage V x is negative, and if V x « -£/ (d) (0)/2 then Ae±i = U {d) (0)(n +1 + n_i)/2 
which corresponds to Fig. 2c of Ref. [14|. That is Eq. (|4.6a| ) completely describes the behavior 



of both minority and majority lines. It makes sense to note that U^(0) at the intermedi- 
ate field is smaller than at high field so the majority line does not shift much, which also 
corresponds to the experiment. 

The last physical phenomena we want to discuss here is the exciton-exciton relaxation 
time. For the estimate of the relaxation time we make use of Eq. ( |5.5| ). The time necessary 
for an exciton to emit or absorb a phonon is around hundreds of ps [f|3| which is much 
longer than short luminescence relaxation times with electrons and hole in the same well 
(a few tens of ps). That means that under the condition of resonant excitation the exciton 
momentum is around that of an exciting photon, K ~ 2 x 10 5 cm -1 . So for concentrations 
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of the order or smaller than 4 x 10 10 cm -2 we can use for an estimate only the second term in 
the numerator of Eq. ( |5.5| ). Since the relevant experiments have been done in multi quantum 
wells structures, where electrons and holes are confined in the same layer, we make use of the 
2D model that we have discussed in the end of Sec. [IV] and estimate the interaction matrix 
element according to Eq. ( |4.10b[ ). In the renormalization factor ( |5.2| ) only the logarithmic 
term can be important and it is ~ 1 - (M /nh 2 )V x ki(Ka) « 8.7 for the exciton radius a « 80 
A. This leads to 1/r ~ 0.12(Hn/ fi) and r ~ 8 ps which is close to the experimental spin 
relaxation time measured by Le Jeune et al. [|12|], Baylac et al. [16], and Amand et al. [|l~7[| . 



Wang et al. |1J| measured the exciton momentum relaxation rate in the concentration 
region between 10 9 to 1.5 x 10 10 cm" 2 in multiple quantum wells of 130 A width. They 
obtained the relaxation rate that grew with the concentration from 0.5 to 2 ps -1 that is 



about an order of magnitude larger than Eq. (5.5) gives. Such a big difference cannot be 
attributed to a deviation for 2D model for exciton. It is likely that a contribution of other 
elastic scattering mechanisms (e.g., surface roughness) was substantial in this experiment. 



VII. CONCLUSION 

We derived the Hamiltonian of exciton gas in quantum wells by the projection of the 
electron-hole plasma Hamiltonian to exciton states and expansion in a small exciton density. 
Matrix elements of the exciton Hamiltonian are expressed in terms of a single exciton wave 
function which is not modified by exciton-exciton interaction and is rather sensitive to the 
geometry of the heterostructure. With the help of the exciton Hamiltonian we estimated the 
blue shift and spin splitting of the exciton luminescence line and their dependence on the 
heterostructure parameters. We also wrote down the Boltzmann equation for excitons and 
estimated the energy and spin relaxation time resulting from the exciton-exciton scattering. 
We succeeded to give an explanation to some recent experimental results that have not been 
explained so far. 
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APPENDIX A: CALCULATION OF MATRIX ELEMENTS 

In this Appendix we present the detailed calculation of the matrix elements that deter- 
mine the exciton-exciton pair interaction. The Appendix contains four subsections where the 

overlap integral A^ ^.^ ^ the direct part of the matrix element Hf^ ^^ ^ 

r («) 

KiSi,K2S2\K 3 s 3 ,,K i s 4 , 



( [3.30| ), the exchange part of the matrix element $ & > (|3.31| ), and the contri- 



bution from the excited states, n^ ex > are calculated. In this calculation we 

omit z-dependent part of the kinetic and potential energy in Ti^eh- This part gives only the 
size quantization energy of free particles which is our energy reference point. 

1. Overlap integral 



The substitution of functions ( |3.24j ) in the definition ( pj.17 ) gives 



Ak iS1i k 2S2; k SS3: k 4S4 — g ^K 1+ K 2t K 3+ K 4 

X \^SleS^S2eS3e^S lh S3 h ^S2hSAhAg 1 ic 2j K 3 ,IC4 ^SleS3 e &S2eSAe$s lh S ih &S2hS3hAK 1 K 2t K 4 ,K-$) ' (^1) 

Sje (sjh) is the electron (hole) spin of the exciton with spin Sj, and 

A KxK2,Kz,Ki = - J dRdr ee df hh exp -%{K X - K^)r ee - i(K l - K 3 )f hh 

i I r> i '"ee ^hh \ , I ft Tee ^hh \ , I fi T"ee Tfih \ , ( ft , T"ee Thh \ / a n\ 

X0^fl+— J— J^^-— J— J^^-— J— J^^fl+— J— J • (A2) 

The overlap integral (|A2|) has an obvious property, Agggg = Ag- g . The re- 
duction of the overlap integral to this form has been made with the help of the relative 
electron-hole coordinates 
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r ee = r el - r e2 , r hh = f hl - f h2 , (A3a) 
R=(f el + f e2 -f hl -f h2 )/2. (A3b) 



The characteristic values of R, r ee , and r^h in integral ( |A^| ) are of the order of the exciton 
radius a. On the other hand, inequality ( |3.2| ) is equivalent to aK <C 1. This means that all 
exponential factors in the integrand can be replaced by unity and 



^K 1 K2,K :j ,K i ~ A = ~ J dRdr^dr^h 



X0 lR+—^—U IR —U IR . (A4) 



The substitution here of the Fourier transform of the wave function, ( [3.28 ), immediately 
leads to Eq. ( ^27|) . 



2. The direct part. 

After the summation over spin variables, Eq. ( |3.30| ) becomes 



xipg s (r e i, r hl )ipg 4 (f e2 , r h2 )dr el df hl df e2 df h2 



2eh 



x^ 4 (^ei, ^1)^3(^2,^2)^61^1^82^*7*2 • (A5) 



The following part of the Hamiltonian ( p. 22 



h 2 v 2 el h 2 v 2 hl _ h 2 v 2 e2 h 2 v 2 h2 

- " + u ^ ~ r «» " " + u ^ ~ ^ (A6) 

is the sum of the Hamiltonians of two free excitons consisting of the pairs (el, hi) and (e2,h2). 
This part of the matrix element gives the sum of the two free excitons energies, Ek 3 + Ek 4 
multiplied by the unit matrix, ( |3.25| ). Here Ek = — e& + h 2 K 2 /2M. For the calculation of 



the other terms of the Hamiltonian, it is convenient to change the in-plane variables of the 
integration to the center of mass coordinate, 
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R c = [m e (f el + f e2 ) + m»(f hl + f h2 )]/2M 



(A7a) 



the distance between the exciton centers of mass, 



p = [m e (f el - r e2 ) + mii(f w - f ft2 )]/M 



and relative coordinates, 



(A7b) 



n = r el - r M , r 2 = r e2 - r h2 



Then the integration with respect to R c results in Eq. (|3.32|) , where 



it-Id + K 2 + K 3 -K± 



X [Ueeflfei - f e2 |) + U hh (\f h i - f h2 \) + M e ft,(|f e i - f h2 \) + U eh (\f e2 - f hl \ 

x 2 (r i ) 2 (r 2 ) dri dr 2 dp 

ii-Kx + Kz-Ks + KA'- 



(A7c) 



(A8) 



+ <5 S1S4 <5 S2S3 / exp 

x K e (|f e i - f e2 |) + Mhft(|r w - r h2 \) + ii e ft(|r e i - r h2 \) + M e/i (|r e2 - r M |)] 
X(j) 2 (ri)0 2 (r 2 )dfidf 2 dp . 

With the help of the momentum conservation that is expressed by Sg 1+ g 2 k 3+ k 4 m Eq. 
( p.32| ), we come up with Eq. ( |3.33|) where 



U {d \q) 



-iqp 



u. 



p + m 



n - r 2 



+Ueh 



P + 



M 

m\\fi + m e f 2 



M 



Uhh 



+ Ueh 



p — m, 
P- 



ri - r 2 



M 

m e T\ + m\\f 2 



M 



X(p 2 (rx)0 2 (r 2 )dfidf 2 dp 



(A9) 



The Fourier transformation with the help of notations, ( |3.35|) - ( |3.36| ), easily reduces this 
expression to Eq. ( |3.34| ). 



3. The exchange part. 

After the summation with respect to spin variables, Eq. ( |3.31| ) becomes 
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K\si ,K 2 S2,K3s :i ,K4S4, 



/ (Pel, rhiWg 2 f e2 , f h2 ) 
x7-(-2ehipK :j (re2, r hl )ijj g 4 (f el , r h2 )df el df hl df e2 dr h2 

/ (f e l ' ^hl ) V£ 2 fe2 , r h2 ) 

*1-(-2eh4>K 4 f e2 , rhi)i>K 3 {rei, r h2 )df e idf hl df e2 dr h2 . 



(A10) 



The fact that these terms correspond to the exchange of two electrons or two holes, so the 
initial and final two excitons are formed from different particles, makes the calculation more 
complicated than the calculation of the direct part. First, we cannot explicitly express prod- 
ucts of the electron and hole spin 5-symbols in the exciton spin as we did in the direct part, 
Eq. (|A5|) . Second, it is harder to separate the single exciton energy Ek from the exciton- 
exciton interaction. Operating by Eq. ( JA6| ) on the wave functions i>gjf e ifh2)' l ^g 3 fe2fhi) 
and il>g 3 (r e i, rh2)4 , K 4 f e2 , ^hi) results in the sum of single exciton energies of the initial state, 
Ek 3 +Ek 4 . However, operating by the same kinetic terms plus different terms of the Coulomb 
interaction, i.e., 



11 +u eh (\f el -f hl \)- 



+ u eh (\f e2 - f h2 \) 



(All) 



2m e 2m\\ ' ' ' 2m e 2m\\ 

on the wave function ?/>*^(r el , r h i)ip*^(r e2 , r h2 ) results in the sum of single exciton energies 
of the final state, E^ 1 + E^ 2 . We write the exchange part in a symmetric form so we add 
and subtract the kinetic terms to the Hamiltonian and reduce Eq. ( |A10|) to 



H 



(*) = 

KlSl,K2S2,K3S3,KiS4 

x E Kl + E K , 2 + E Ks + E Ki + 



^S le S4e^S2eS3 e ^S lh S3 h Ss 2 hS4h j ^ K\ 

{f\fh\)^\ 2 f e2 ,f h2 ) 



^li , ^ 2 vL . n 2 v 2 e2 , n 2 w 2 h2 



2m P 



+ 



2mw 



+ 



2m P 



+ 



2m\\ 



+Uee(\r e i - r e2 \) + u hh (\f hl - r h2 \) 
x^ 3 (r e2 , T h x)ij}gJf A , r h2 )dr el dr hl dr e2 df h2 



^S lh S4 h ^S 2 hS3h / fel, rhlW^ fe2, r h2 ) 



X 



Er x + E K2 + E Ka + E Ka + 



*M , h 2 V 2 hl , h 2 V 2 e2 , h 2 V 2 h2 



2m P 



+ 



2m\\ 



+ 



2m P 



+ 



2m\\ 



+u ee (\f el - r e2 \) + u hh {\f hl - r h2 \) 
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xipg (f e i, r h2 )if;g (f e2 , r hl )dr el df hl dr e2 df 



h2 



(A12) 



Convenient variables for the calculation of the integral are coordinates R, r ee and fhh, Eq. 
), and the center of mass of the system R c , Eq. ( |A7a| ). After integration with respect to 



R c we obtain 



H 



(x) 

KiSl,K 2 S 2 ,K 3 S 3 

X 



,i? 4S4 - g 5 K 1 +K 2 ,K 3 +K i {^Kl + E «a + E K S + E K 4 



h 2 K 2S 

AM 



(A13) 



1 



X 



where 



#1^2,^3,^4 



dRdf ee df hh 



x exp 



-i{K x -K 2 ) 



Tee + r hh 
2 

ft 2 



Tee ~ T hh 



h 2 h 2 

Uee{ree) + U hh {r hh ) + '—V 2 + — V* e + — V\ h 

4/z m e 



x exp 



-i(K 3 - K 4 ) 



Tee ~ r hh 



R 



Tee + r hh 



R 



Tee ~ r hh 



R + 



Tee + r hh 



(A14) 



From Eq. ( |AT^ ) it follows that 



Tj{ x ) _ Tj( x ) 

K!K 2 ,K 3 ,K 4 K 3 ,K 4 ,KiK 2 



(A15) 



The matrix element, U;~ \-i r - r , , is simplified with the help of the low exciton energy 

' K 1 K 2 ,K 3 ,K 4 ) r r &j 



assumption, (|3.2j ). This is done in a similar way to the simplification of A^g^^g^, in 



Appendix Al. The result does not depend on the excitons momenta, 



U^- - - = = — I dRdr druu 



x0 R + 



r h h 



R- 

" h 2 



Tee ~ T hh 



h 2 h 2 

Uee(r ee ) + u hh (r hh ) + —V 2 H V 2 e H V 2 hh 

Ajjt m e m h 



X(f>[R 



Tee + r hh 



R + 



Tee + r hh 



(A16) 
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This integral takes a more compact form if it is expressed in Fourier transform of the wave 
function, 

^ 2 f 4.4 2 dq f 2 , 2 /\- /\- - im dq\dq2 



u(x) = Jjrt q2 jSy-j « [u - (|(fl " ^ ] + Uhh (l * ~ m W ' (A17) 



This last integral can also be written in another form, 



U(x) = jjrt f J§yt- lKe(r)+u hh (r)} 



9 (2tt) 2 



2 



df . (A18) 



Integral ( |A16| ) is of the order of h 2 / \i. Since e b ~ h 2 / /ia 2 , the second term in Eq. (|A13|) is 
of the order of U^ x > /S ~ £ b (cL 2 / S). According to Eq. (|A"4l) , A ~ a 2 , hence the first term in Eq. 
( |A13 ) is of the same order as the second. The exciton kinetic energy in the first term can be 



neglected because its contribution to H^) - r , r , is of the order e b (a 2 /S)(h 2 K 2 /M), 
i.e., contains the product of two small parameters. As a result Eq. (|A13|) is reduced to Eq. 



4. Contribution of excited states 



For two excitons the definition, (|3.14 ), can be written as 



H t::L,,K^= E (KiS 1 K 2 s 2 \H + 2e b \q 1 a 1 s' 1 q 2 a 2 s' 2 ) (A19) 



<3l.<32. a l' a 2 
s l' s 2 



/ ^2 2 ^2 2 \- x 

x f E ai + E a2 + —± + — j- + 2e b \ (qiaxs'^a^lH + 2e b \K 3 s 3 K 4 s 4 ) . 

The matrix elements of the Hamiltonian between the ground state and excited states can 
be separated into the direct and exchange parts [compare Eqs. (|3.29|) - (|3.31|) 1, 



(K lSl K 2 s 2 \H + 2e b \q 1 a 1 s[q 2 a 2 s' 2 ) = H^^^R^, K 2 s 2 ) + H^^^^K^s*, K x a x ) 

-HgL^t&si, K 2 s 2 ) - H^ 2a2S ,(K 2 s 2 , K x s x ) . (A20) 

In the calculation of these matrix elements it is convenient to separate the sum of two 
single-exciton parts of the Hamiltonian, ( |A6| ) , and to operate it on the left bracket [similar 
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to the calculation of ~ - -. and H y i' -. - - 1. Then, after the neglect of 

KiSl,K2S2,K3S3,K4,S4 KiSi,K2S2,K3S3,K4S4 J u 

the exciton kinetic energy, 



r(x) 



I df el df e2 df hl df h2 

(J e i(J e 20- hl CT h 2 

x *^i,.i(^i» r hl , a hl )tyg 2 S2 (r e 2, a e2 ; f h2 , a h2 )U(f el ,f e2 , f hl , r h2 ) 
x^olais', (^ei, o- e i; r hl , cr hl )^ff 2a2S ' (f e2 , cx e2 ; r h2 , cr fo2 ) , 



(A21a) 



where 



E 



(T el (T e 2(T, ll O-( l 2 



df el df e2 df hl df h2 



X^f*g i si (f el , a el ] f hl , 0-w)*^ 2jS2 (^e2, 0" e2 ; r h 2, O h2 )V{f eX , f e2 , f hl , f h2 ) 

><^g iaiS [(r el ,a el ]r h2 ,a h2 )^^ a2S ' 2 (r e2 ,a e2 -,r hl ,a hl ) , (A21b) 



U(r el ,r e2 ,r hl ,r h2 ) = u ee (\r el - f e2 |) + ww^lf/a - f ft2 |) 

+ Uehd^el - ^hal) + U eh (\f e2 - f hl \) 



(A22) 



After the substitution of the exciton wave function, Eqs. fl3l]) - flU), in Eq. ( |A2l| ), the 
spin factors are easily separated, 



H 



qiais 1 q2Ct2S 2 
r(») 



/(A'iSi, K 2 S 2 ) = 5 S1)S ; ■ +s; '5 S2 V .+.(.#. 



(d) 

e2~ r!s h2~~q.\ a \m a 1 



S el"T' , h2 * 2 ' i e2^ A hl <?l a l> <?2"2 



(i?!,K 2 ) 



(A23a) 
(A23b) 



and 



(d) 

q\ct\m<±2 



(K U K 2 ) = — J df el df e2 df M df h2 U(f el ,f e2 ,f hl ,f h2 ) 
m e r el + m\f h i _ rrief e2 + m\\f h2 



H, 



(x) 

gioi,g202 



x exp 
x exp 

x exp 
x exp 



M M 
m e f el + m\\f hl m e r e2 + m\\r h2 



(A24a) 
<t>{\r e i -r h i |) 0(| f e2 -f h2 1) 



M 



+ «92" 



M 



0oi(^*ei - r hl )<f) a2 (r e2 - f h2 ) 



S 2 



df el df e2 df hl df h2 U (f e i, f e2 , r hl ,r h2 ) 



. r > m e r e i + miir/ji . - m e r e2 + m\\r h2 

-iK\ — — - iK 2 — — 

M M 

m e r el + m\\r h2 m e r e2 + m\\f hl 



(A24b) 
4>{\f el - r hl \)(p(\f e2 -r h2 \) 



M 



+ iq2- 



M 



0oi(^ei - r h2 )(p a2 (r e2 - r hl ) . 
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The integrals are simplified with the help of new integration variables. For the direct part 
they are given by Eq. ( |A7| ) and for the exchange part they are the two-exciton center of 
mass, R c , Eq. ( |A7a| ), and the relative coordinates, Eq. (|A3|). After the integration with 
respect to R c 

Hfl lAa2 {Ku K 2 ) = i <fc 1+ * aA+a Dg^Ku K 2 , q u q 2 ) , (A25a) 
H^ 2a2 (Ki,K 2 ) = I 6g 1+ g aA+ ^ D^itu . (A25b) 



In the remaining integrals, D$ and , it is possible to neglect K x and K 2 which 

due to 5-symbols in Eq. ( |A25|) means also that q 2 = —q\- This allows us to simplify the 
notations, D$ a2 (K u K 2 , q 1: q 2 ) = £>gU<Zi), D^ a2 (K u K 2 ,q x ,q 2 ) = D^fa). For further 
simplification it is convenient to substitute the Fourier transform of the exciton functions, 

<f>a4= J e-^ct> a {f)dr . (A26) 

The calculations are a bit cumbersome but straightforward and they lead to Eqs. ( |3.40| ) - 

ra>- 



APPENDIX B: SPIN SUMS IN THE EXCHANGE TERM 

In this appendix we reduce the spin sum in the exchange term of the many exciton 
Hamiltonian, Eq. ( [3.52j ), to a simpler form. According to the definition, Sj = s e j + Shj, and 
there is one to one correspondence between the electron and hole spins on one side and the 
exciton spin on the other. Therefore 

Sl,S2,S3,S4 

XC k,S2 C k,Sl C ^l-9. s 3 C ^2+g,S4- 



/y {$Sel,Se4,$s e 2,Se3$s h i,S h3 fis h2 ,S h 4 + ^S e l,S eS ^S e2 ,S e 4^S hl ,S h4 ^S h2 ,S 

e4,Shl, s h2, s h3, s h4 

(y— c ~* c ~* 

'K 2 ,s e2 +s h2 K 1 ,s el +s hl K!-q,Se3+s h3 K 2 +q,s e 4+s M 
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Sel,S e2 ,S e 3,S e 4,S hl ,S h2 ,S h3y S h 4 

C t_ — (Bl) 



^ K 2 ,s e2 +s h2 K lt s el +s hl V K 1 -q,s e2 +s hl K 2 +q,s e i+s h2 Kt-q,s e i+s h2 K 2 +q,s e2 +s h 

s el, s hl, s e 2 ,Sh 2 
= S} + S 2 + S3 + S4 , 

where 

S i = 2 E c /? 2iSe+s , i c /? liSe+Sh c i? 1 -g, Se + Sh c /?2+g,^+ Sh > ( B2a ) 

^ 2 = E C i? 2 . Se+Sh C i? 1 _ Se+Sh ( C i?i-g -s e +s h C K 2 +q,s e +s h + C £ 1 -£ Se+Sfc C £ a+ £- ae+Sh ) ' (B2b) 
S 3 = E C £ 2lSe + Sh C /?i,s e -s h ( C i?i-<f,s e +s fe C i?2+?>e-Sfc + C Ki-q,s e -s h C K 2 +q,s e +s h ) ' (B2c) 
^ 4 = E 4? 2 .s e + Sh C /?i -« e -S h ( C i?l-<?,-S e +Sh C i?2+g,S e -S h + ^i-^Se-SA^S+g-Se+Sfc) ' (B2d) 

Here we use the fact that the both the spin of the hole and the spin of the electron have 
only two values, s e = ±1/2, Sh = ±3/2. Hence, the only possible values of are ±s^i 
and of s e2 are ±s e i- Sj describe all the possible combinations that appear in Eq. (|B1[). 



The next step is to replace the quantum numbers of the electron and hole spin with the 
quantum numbers of the excitons spin. Here we use the one to one correspondence between 
the electron and hole spins and the exciton spin 

5 1 = 2 E C k,s C k,s C Ki-€,s C K2+q,s » ( B3a ) 

s 

5 2 + S3 = ^ C 7? 2 , S2 C i? llS1 { C Ki-q,S! C K 2 +q,s 2 + ^-(f^ifa+Jsi) 

Sl,S 2 

~ 2 E C ^ 2 . s C i? liS C /?i-g, S C i?2+(?,s ~ E C k,s C k -a { C Ki-q,s C K 2 +q,s + C K 1 -q,s C K 2 +q. -s) ' (B3b) 

S 4 = E ct Clt Cff-. * „„C 



E4?2,s C k -s ( c i?i-g-s c i? 2 +g> + c i? 1 -<f- s c i? 2 +<f,s) " ( B3c ) 



The exchange term in the many exciton Hamiltonian, Eq. Q3.52j ), has now the following 
form 

H {X) = ^ E (S 1 ±S 2 + S 3 ±S 4 ), (B4) 

Ki,K 2 ,q 

Now we replace the summation over q with the summation over Q = K\ — K2 — q in the 
second term in both parentheses in S 2 + S3 and see that it becomes equivalent to the 
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first term in the same parentheses. The same can be done with the two terms inside the 
parentheses in S 4 . So the exchange term is 



it(x) = Y^. v 

K u K 2 ,q 



2E 

Sl,S 2 



C k,s 2 C k,si C ^i-9>i C ^2+g>2 



(B5) 



^^ C ^K 2 - s C ^K 1 , s C K 1 -q,s C K 2 +q-s+ C K 2 ,- sl C K ljS1 C iCi-q~,S2 C iC2+q, 



S2 



Sl,S 2 



APPENDIX C: FOURIER COMPONENTS OF THE INTERACTION 
POTENTIAL AND THE WAVE FUNCTION IN QUANTUM WELLS 



In this appendix we present some formulae necessary for the calculation of Hamiltonian 
matrix elements. 

The expressions for Uij(r) that are obtained with the help of Eqs. (|2.6|) and ( |4.5| ) are 
quite complicated. But their Fourier components have a relatively simple form, 

2vre 2 (8vr 2 ) 2 sinh (qL e /2) sinh (qL h /2) 



Ueh(q) 

Uhh(q) 

U ee (q) 



qK L h L e g 2 (47r 2 + L 2 g 2 )(4vr 2 + L 2 g 2 ^ 



-q(w+L h /2+L e /2) 



2ne 2 
qK 
2vre 2 



+ 



qU. 



qL h ' 4vr 2 + g 2 L 2 



qK 



+ 



qL t 



32tt 4 (1 - exp {-qLh)) 

q 2 L\ (4tt 2 + g 2 L 2 ) 2 
32tt 4 (1 -exp(-gL e )) 



qL e An 2 + q 2 L 2 e 



(Cla) 
(Clb) 
(Clc) 



g 2 L 2 (4yr 2 + q 2 L 2 e ) 2 

Each of these functions is singular at q = 0. But the direct interaction potential which 
contains the sum of them is finite, 



2%e 2 r (2 5 . , 



K 



(C2) 



This expression leads to Eq. 

The Fourier transform of single-exciton wave function, ( }4.6|) , necessary for the calculation 
of UW(q) and V x , is 

-3/2 



b(b + r ) 2b \4b 2 q 



1 + r °V46 2 + ^ ' 6XP 



~ r lW 2+q2 



(C3) 
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FIGURES 

FIG. 1. The dependence of the exchange coefficient, V x = — 2\eb\A, on the separation 
between the wells, d = w + L. The continuous line is for fixed barrier width, w = 42 A, and 
changing wells width. The dashed line is for fixed wells widths, L = 70 A, and changing the barrier 
width. 

FIG. 2. The dependence of the exchange coefficient, V x , on the separation between the wells, 
d = w + L, where w = 10 A and L is changing. V x becomes positive for a small enough separation 
between the wells. 

FIG. 3. The dependence of the coefficient, Vj,, on the separation between the wells, d = w + L. 
The continuous line is for fixed barrier width, w = 42 A, and changing wells width. The dashed 
line is for fixed wells widths, L = 70 A, and changing the barrier width. 

FIG. 4. The dependence of the energy splitting parameter, V es , on the separation between the 
wells, d = w + L, where L = 70 A and w is changing. 
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